We reconstruct a 3+1 formalism of general relativistic electromagnetism, and derive 
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> 1 ABSTRACT 

O . 

in ■ 
\o ■ 

the equations of motion of charged particles in the pulsar magnetosphere, taking account 
of the inclination between the rotation axis and the magnetic axis. Apart from the 
■ previous works where space charge is evaluated by assuming the flow velocity being the 

light speed, we analyze particle motion in the polar cap, finding that gravity changes 
q . significantly its dynamics and the condition for acceleration. 

a: 

Subject headings: acceleration of particles, gravitation, magnetic fields, pulsars 

> : 

s_j ■ 1. Introduction 

d; 

The origin of radio emission from pulsars remains a great mystery. One of the likely sce- 
narios is that particles are accelerated along open magnetic field lines and emit 7-rays which 
subsequently convert into electron-positron pairs under a strong magnetic field. A combination 
of the primary beam and pair plasma provide the radio emission mechanism (Melrose 2000). The 
process of particle acceleration has been intensively studied in flat spacetime (Fawley et al. 1977; 
Arons & Scharlemann 1979; Mestel 1981; Shibata 1991; Shibata 1995; Shibata 1997). The field- 
aligned electric field is driven by deviation of the space charge from Goldreich-Julian charge density 
(Goldreich & Julian 1969), which is determined by the magnetic field geometry. Therefore, gen- 
eral relativistic (GR) effects on the field geometry are crucial for the formation of the field-aligned 
electric field and particle acceleration. 



Muslimov & Tsygan (1992) initiated the GR analysis of electromagnetic fields around a pulsar; 
they solved Maxwell equations on the assumption that the particles move with the light speed. 
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Shibata showed, however, that their assumption of constancy of particle speed was not always true 
(Shibata 1991; Shibata 1997): the particle motion and space charge show oscillatory behavior even 
on field lines curving toward the rotation axis if the assumed current density is sub-critical. Mestel 

(1996) extended Shibata's equations of motion to include GR effects, though his analysis is not 
complete. Recently Harding and Muslimov (2001) extended Muslimov-Tsygan's analysis to include 
the pair-creation dynamics. 

In this paper we derive the basic equations more rigorously and generally: we depart Muslimov- 
Tsygan (and Harding-Muslimov) model by solving equations of motion for the particles, and correct 
Mestel's equation and extend it to include the inclination angle between the magnetic axis and the 
rotation angle. Next we solve the electric field together with the particle motion for the region near 
the magnetic pole just above the surface. It turns out that the effect of gravity on particle dynamics 
cannot be ignored, and hence the condition for acceleration is also modified. As pointed by Arons 

(1997) , local deficit of space charge takes place on most of the field lines. In particular, particle 
dynamics on the "away" field lines from the rotation axis is changed qualitatively by gravity if the 
current density is less than the critical value. It is also true that oscillatory solution exists even 
with GR effect. We find that some oscillatory solutions disappear above a certain height. 

This paper is organized as follows. In section 2, we reconstruct a 3+1 formalism of GR electro- 
magnetism. In section 3, we derive equations of motion of charged particles around a magnetized 
rotating star. In section 4, solving the equations of motion, we analyze particle acceleration in the 
polar cap. Conclusions and discussions are given in section 5. We use the unit c = 1. 



2. 3+1 Electromagnetic Equations: General Formalism 

A 3+1 split of Maxwell equations were made by Landau <fe Lifshitz (1975) and Thorne & 
Macdonald (1982). Landau & Lifshitz derived 3+1 Maxwell equations in a straightforward manner, 
which we essentially follow. They introduced a metric 1 , 

ds 2 = g^dx^dx" = —h(dt — gidx 1 ) 2 + r )ijdx l dx^ (1) 

h = -goo, 9i = goi/h, i j = g ij , (2) 

which is formally in a 3+1 form. However, it is hard to understand how the entire 4D spacetime 
{x^-,9iiv) is decomposed into 3D spatial hypersurfaces (x l ,qij) because each hypersurface does not 
correspond to t =constant. 

On the other hand, Thorne & Macdonald adopted the metric used in the ADM formalism: 

ds 2 = g^dafda? = -a 2 dt 2 + q ij {dx i + f3 i dt){dx j + pdt), (3) 



1 In this paper Greek letters run from to 3, and Latin letters run from 1 to 3, contrary to the convention of 
Landau & Lifshitz (1975). 
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a 2 = -1/500, A = 90i, Qij = 9ij, (4) 

where a and (3 are called the lapse function and the shift vector, respectively. With this metric one 
can understand the concept of the 3+1 decomposition easily, as we describe below. However, they 
adopted the "congruence" approach and wrote down the equation with coordinate-free language 
such as expansion and shear of the fiducial observers. Although their equations are powerful for 
some theoretical arguments, one cannot avoid using coordinates to solve them in most cases. 

For the purpose of our study, we combine and modify the two methods: we adopt the ADM 
formalism to decompose a spacetime and electromagnetic fields into a 3+1 form, and write down 
3+1 electromagnetic equations in a straightforward manner. 

First, with the metric (3), we introduce the fiducial observers with the 4- velocity, 

«„ = (-<*, 0,0,0), (5) 

and the spacelike hypersurface £ which is orthogonal to at each time t. Because £(£) is charac- 
terized by a t = constant hypersurface, x % and qij are the 3-coordinates and the 3-metric of 
respectively. We define the projection tensor, which projects any 4-vector or tensor into £(i), as 

= gf" + . (6) 



Electromagnetic quantities are originally in the 3+1 form, i.e., 3-vectors and scalars: electric 
field E, magnetic field B, charge density p, current density J, scalar potential 4>, and vector 
potential A. We define those quantities on £(i), i.e., those quantities are supposed to be measured 
by the observers with the velocity u^. Note that the "ordinary" components of a vector P are 
neither the covariant components Pi nor the contravariant components P l (Weinberg 1972). If the 
3-coordinates x l are orthogonal, i.e., = h^Sij, the ordinary components P^ are 

P {i) = h {l) P = P l /h {i) . (7) 

As long as we adopt the 4- velocity (5), which specifies S, we can convert any 3- vector P on S is 
converted into a 4-vector P M in a natural way, 

= ( 0) pi). (8) 

It is easy to confirm that P M is the vector on S: P^u^ = and P^h 1 ^ = P 1 . 

Let us reconstruct the electromagnetic field T^ v ', the charge-current vector J 11 and the 4- 
potential from the quantities defined on E and (Thorne & Macdonald 1982): 



(9) 
(10) 
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where e^ vXa is the Levi-Civita tensor with e 0123 = 1/y/^g. We can invert these relations as 

& = F v u v , = - X -^u v T Xc , (11) 

P = - J* = KJ V , (12) 

cj) = -4%, A = ti v A v . (13) 

The 4D expressions of Maxwell equations and the continuity equation are 

= P { ^ ]X] = 0, J^ = 0. (14) 

Substituting (9) and (10) into (14), we obtain 3+1 Maxwell equations: 

divB = 0, (15) 

-l-dtiy/qB*) + [vot(aE + (3 x B)\ = 0, (16) 

div£ = 4vrp, (17) 

— ^(v^*) + [rot(aB + /3 x = 47raJ*, (18) 
v ^ 

4=^(V9P)+dM«^) = 0, (19) 
where cross product, divergence and rotation are defined as 

[Q x Py = e ijk QjP k , d\wP = ^-di{^P l ), [rotP] 1 = e ijh djP k , (20) 

and e yfe is the spatial Levi-Civita tensor with e 123 = l/y/q. 

To complete 3+1 electromagnetic equations, we rewrite the equation of motion of a charged 
particle. We denotes its mass by m, its charge by e, its proper time by r and its 4-velocity by 
= dx^/dr. Then the equation of motion is written as 



m -7- +T>V (21) 



(' lrl ' . ,-,< 

where T is the affine connection. We define a scalar 7, the 3-velocity vector V 1 and its norm V as 

7 =- v ^ u -av°, V i = —^ = - +PV , V^q^VK (22) 

7 7 

From the normalization v^v^ = — 1 and (22), we find 7 = — V 2 , which implies that the scalar 

7 is nothing but the Lorentz factor. With these quantities and the 3+1 electromagnetic quantities 
we rewrite the the equation of motion (21) as 

A. {7 (_ a + p . v )} _ f. ( doa + ^V^d m + V • do/3) = • V, (23) 

dr a \ 2 / ma 

rf /..t^ 7 2 / Q . , a T m,in 1 t 7* a ^\ e 7 



(7^) - - (-dia + ^V*^-* + V ■ diP) = ^- {aE + V x B - (E ■ V)(3} 1 . (24) 
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3. Electromagnetic Equations in Pulsar Magnetosphere 

We assume that the spacetime outside a neutron star is stationary and axisymmetric. Because 
gravity of charged particles is negligible, the outer spacetime is described by Kerr metric. We adopt 
the spherical coordinates (t, f, 9, up) so that 9 = accords with the rotation axis. In the slow- 
rotation where we take up to the first order of the angular momentum of the star L, the spacetime 
is described by the metric, 

f 2L 
a 2 = F(r) = 1 ~j, = = fa = diag(-F _1 (f), f 2 , f 2 sin 2 #), (25) 

where r g is the gravitational radius. 

We consider the two-step coordinate transformations. First, we introduce the "corotating 
coordinates" as 

tp — ><p- nt, (26) 

where Q, is the angular velocity of the star (Landau & Lifshitz 1975). The shift vector is, accord- 
ingly, transformed as 

P? — ► j& = Q - to. (27) 

Second, we move on to the "magnetic coordinates" (t, r, 9, (p), where the magnetic dipole axis 
accords with 9 = 0. Supposing the magnetic axis is at 9 = x, <p = 0, the coordinate transformation 
from the corotating coordinate to the (corotating) magnetic coordinate is given by 

x = f sin 6* cosy?, y = r sin 6* sin z = rcos9, 

x = xcosx + -Ssinx, y = y, z = — xsinx + icosx, (28) 
x = r sin 9 cos 99, y = r sin 9 sin ip, z = rcos9. 

The shift vector in (25) is, accordingly, transformed into 

(3 r = 0, (3 d = -sm X smp(n- uj), = (cosx - cot6»cosv9sinx)(^ - u), (29) 

while the lapse function and the spatial metric remain unchanged: 

a 2 {r)=F(r), qij = diag(F _1 (r), r 2 , r 2 sm9). (30) 



In this (and any) axisymmetric spacetime the observers with the velocity (5) is called "zero- 
angular-momentum observers" (ZAMOs) because of vanishing their angular momentum per unit 
mass: u^d/dipY = (Thorne & Macdonald 1982). 

Now we analyze the basic equations (15)- (19), (23) and (24) with the metric (29) and (30). 
First, we suppose that the magnetic field distortion due to the external currents is negligibly small. 
Then the magnetic field is governed by (15) and (18) with J = 0: 



d\vB = 0, iot(aB) = 0. 



(31) 
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We adopt the dipole-like solution of (31) (Ginzburg & Ozenoi 1965): 



B (r) = 



B {9 ) = 



fir) = 



B ( 



1 2 B = ^^ T [-2f(r) + 3F- 1 (r)] 



fir*) 



2 fin) 



sin 9 



COS I 



(32) 



3r f _ 



3 



sin I 



3r„ 

1 + 



(33) 
(34) 



where r = denotes the star surface, and the approximate expressions are obtained by expanding 
in powers of r g /r or r g /r* (Konno & Kojima 2000). 



Equation (16) is satisfied if there exists a potential <J> such as 

aE + (3 x B = -V$. 
Substituting (35) into the other of Maxwell equations (17), we obtain 



div = 4tt(p - p GJ ), 4irp GJ = -div — x B , 



a J 



a 



where pcj is the Goldreich-Julian density. Each term is written explicitly as 



div 



\ a 
Airpdj 



F d 
dr 



dr 



+ 



1 



^9d9 ( Sme M 



+ 



1 d 2 $ 



-aB ■ rot— ~ 

or 



B 



(r) 



v^Fsintf 



d(sin 2 00*) dp 



39 



dtp 



sin 2 6 dip 2 
Fsm9 d (r 2 pv 



(35) 



(36) 



(37) 



(38) 



Let us consider the flow of charged particles. In the case under consideration, magnetic field 
dominates over electric field, where particles are at the lowest Landau level and the inertial drift 
motion across the magnetic field is negligible. We may, therefore, regard that particles go along 
with magnetic field lines: 

V = kB, J = P V = pkB, (39) 
where n = V / B = \ V\/\B\ is a scalar function. 

Then continuity equation (19) implies (Muslimov & Tsygan 1992; Mestel 1996) 

= div(aJ) = di-v(apKB) = B ■ V(apn), (40) 

which reads 



apK 



FJ 
B~ 



constant on magnetic field lines. 



(41) 
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Equations (23) and (24) give the motion of charged particles. If we assume B v = = 0, 
equation (35) and (39) give V<p = = 0. Using the relation 

E . y _ _Zi . y _ _ _I (1^ + 1) |£ _ 1 *, (42) 

a a a \T dr a J dx l ^/FT dr K ' 



we rewrite one of the equations of motion (23) to 



-(VFT) = ^ a * (43) 
At mF At 

The other equation (24) need not be solved because the spatial trajectory is determined by (39). 

Finally, let us consider boundary conditions. We assume that the star crust is a perfect 
conductor, and hence particles on the surface do not suffer Lorentz force. Equations (24), (35), 
(39) and (42) reduces this ideal-MHD condition to V<J> = 0. We also assume that there exist closed 
magnetic lines, where the ideal-MHD condition holds. Thus open field lines in the polar region have 
another boundary, defined by "last open field lines", 9 = 9 c (r). Therefore, our boundary conditions 
are <E> = V<1> = on r = r* and on 9 = 9 c (r). 



4. Particle Acceleration in the Polar Cap 

We are interested in the "polar cap" region, 9 < 9 c (r), where particle acceleration may occur. 
For simplicity, we restrict ourselves to the region near the magnetic pole just above the surface, 
where the following approximations hold: (i) 9 c (r) ~ constant, (ii) 9 <C 1 (we take its first order), 
and (iii) d/dr along the particle trajectory « d/dr. Therefore, we can apply our equation below 
only to the region elongated less than the polar cap radius: r — r* < r9 c . 

As for the non-radial direction, we expand as 

<S> = Y,®{r)Y lm {9,v)- (44) 

l,m 

We only take the mode of the polar cap scale, I « tt/9 c . 
The above assumptions simplify (36) as 

~ — Tr [ r V ) + 7^ = \V ~ PGJ ) ' (45) 

4vrp GJ = { ( X ~ ^) cos X + 9Q cos </? sin x} , (46) 

Q _ r 9 3" . r g u 3 / Sr g a; \ 3 llr g _ 3w 

V 2r 2fi rfi 2fF \ 2r 2Q J ~ 2 r 40' 1 ' 
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where B = \B\ = B^ + O{0 2 ) 2 . Equation (43) describes the time-evolution of particles; however, 
if particles make a stationary flow, we can read their spatial configuration by replacing the time r 
with a spatial valuable r: 

d , i— x e 1 d$ . . 

dr mF dr 

Using the normalized variables (Shibata 1997), 



2kJF{s)J{s) if ^ k , , j2QB e 

= nnf \ = constant, 4> {s ) = (s , a = */ £-r, 50 

l!i)(s) m V mr 



j(s) = v . > = 1 — cosx + 6>(s)Q(s)cos(^sinx, (51) 



Q{s)B(s) V n 

we rewrite (45) and (49) as 



1<^> ^ f£ <»> 

The particle trajectory is given by integrating d6/dr = B r /B d , resulting in 

9(s) « 0*\f^- for < 1, (54) 

v ^* 

where * denotes a value at the star surface r = r* for any variable. In the following analysis we 
fix some of the parameters: r* = 10km, r g = 2km, VL = 27r/0.1s, uj*/VL = 0.1, 9 C = ^JVtr^jc « 
0.046 (vr/fl c «= 69), X = 30° and 7* = 1.00001. 

We can argue equation (52) on the analogy of particle motion by regarding <p as position and 
s as time. At a given point (or time) j > j(s), the "force term" (LHS) of (52) always positive; <f> 
and 7 increase acceleratingly. If j < j(s), on the other hand, the force term takes both positive 
and negative, depending on V, and hence oscillating behavior is expected. If we set r g = uj = 
(no gravity) and x = (no inclination), we reproduce Shibata's result, j = 1 (Shibata 1997). 
The expression of j in (52) indicates that both effects of gravity and of inclination reduce j, or 
equivalently, the critical value of j. 

Numerical solutions of (52) and (53) in Figs. 1-3 verify the above arguments. It is apparent 
that gravity reduces the critical value of j: monotonic increase of 7 appears with lower values of j 
that in non-GR case. 



2 Equation (47) is slightly different from the expression obtained by Muslimov and Tsygan (1992) and Harding 
and Muslimov (2001): 

W 2r 2tt 2fF \ 2r 2Q. J y ' 

We believe that there is an error in their calculation. Nevertheless this difference does not change results practically; 
because we have assumed 6 <C 1, all terms but 2/3 (the dominant one) in Q are negligible in equation (46). 
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Let us consider the dynamics for j ps j* carefully. To see the behavior analytically, we expand 
j with £ = s/s* — 1 up to the first order, resulting in 



This implies that j increases or decreases according to the sign of C. For the particles on the 
magnetic axis (9 = 0), j is an increasing function due to the term in GR case while it is 

constant in non-GR case. This gravitational effect is illustrated in Fig. 1: contrary to non-GR case 
(a), in GR case (b) oscillation for j < j* does not continue to infinity but the particle velocity 
becomes V = at some point. Although the dynamics is independent of 7* in non-GR case, the 
ending point V = (7 = 1) depends on 7* in GR case: it becomes further for larger 7*. 

Dynamics of particles off the magnetic axis (9 7^ 0) is more interesting. In the absence of 
gravity (oj* = 0), the sign of C depends simply on whether particles move "toward" to (cos (p > 0) 
or "away" from (cos ip < 0) the rotation axis, as claimed by Shibata (1997). He showed, for example, 
that particles are accelerated after oscillation along "away" field lines if j is slightly less than 1. This 
behavior is reproduced in Fig. 2(a). Gravity changes this behavior drastically. If we take uj*/VL = 0.1 
typically, C becomes positive even on "away" field lines unless 6* > 0.4| sec <p\ cot \ > 0.4 cot 
Unless x is large enough, in most small-0 region, j simply increases and acceleration after oscillation 
cannot occur. This argument is confirmed by the numerical result in Fig. 2(b). 

For the particles moving "toward" the rotation axis, cos <p is positive, and therefore j always 
increases along the field lines. The results in Fig. 3 are also understandable. Among the solutions 
in Fig. 3(b), the solutions with j = 1.003j* could be realistic because particle energy becomes so 
large that electron-positron pairs are created with finite electric field. 



We have reconstructed a 3+1 formalism of general relativistic electromagnetism, and derived 
the equations of motion of charged particles in pulsar magnetosphere. Our basic equations are 
the correct and generalized version of those of Muslimov & Tsygan (1992) and of Mestel (1996) 
in the sense that the equations include arbitrary current density, arbitrary velocity which satisfies 
the equation of motion, and arbitrary inclination angle between the rotation axis and the magnetic 
axis. 

We have solved our equations of motion together with the electric field structure along the 
magnetic field for the region near the magnetic pole just above the surface in our approximate 
method. We have found that the effect of gravity on particle dynamics cannot be ignored and hence 
the condition for acceleration is also modified. In particular, particle dynamics on the "away" field 
lines is changed qualitatively by gravity if the current density is sub-critical. 

Monotonic increase of 7 with the boundary condition E ■ B = at the stellar surface and at 




(55) 



5. 



Conclusions and Discussions 
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infinity is found with GR effect, and the achieved Lorentz factor is significantly larger than the 
case of non-GR (Muslimov & Tsygan 1992). With sub-critical current densities, we find oscillatory 
solutions even with GR effect. Since we take the current densities as a free parameter, the local 
accelerator model can be linked with the global model by adjusting the current density. 

We have also found that some oscillatory solutions disappear above a certain height. This 
may indicate an intermittent flow of charged particles rather than a steady outflow if the imposed 
current density is smaller than a threshold. This may cause pulse nulling. 

As pointed by Arons (1997), an important GR effect is local increase of j, which appears on 
most of field lines, no matter whether large scale field curvature of magnetic field lines is "away" or 
"toward". He suggested that this fact can be responsible for the axisymmetric radio emission about 
the magnetic axis. If the current density is taken as a free parameter, the super-critical current 
density in the polar annuli can also be responsible for the distribution of radio emitting region. We 
shall extend our analysis to a more global region to construct a pulsar model in the subsequent 
work. 

N.S. thanks Yasufumi Kojima for discussions, and also Yukawa Institute for Theoretical 
Physics, where a part of this work has been done as a visitor program. This work was partially 
supported by the Grant-in-aid for Scientific Research Fund of Monbukagakusho (Nos. 13740139 
and 12640229). 
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Fig. 1. — Spatial distribution of 7 along magnetic field lines. We choose 6 = 0. (a) represents 
solutions for the flat spacetime with r g = to = 0. 7 oscillates or increases monotonically according 
to j. (b) represents GR solutions. In the case of j < j* oscillation eventually stops with V 
approaching zero, in contrast with the non-GR case. 
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Fig. 2. — Spatial distribution of 7 along "away" field lines. We choose 6** = 5°, ip = 180°. (a) and 
(b) represent flat and GR solutions, respectively. Gravity changes particle dynamics for j < j*: 
acceleration after oscillation does not occur. 
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Fig. 3. — Spatial distribution of 7 along "toward" field lines. We choose 6* = 5°, 9? = 0°. (a) and 
(b) represent flat and GR solutions, respectively. 



